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Neutron stars are the most dense objects in the observable Universe and conventionally one 
uses nuclear theory to obtain the equation of state (EOS) of dense hadronic matter and the global 
properties of these stars. In this work, we review various aspects of nuclear matter within an effective 
Chiral model and interlink fundamental quantities both from nuclear saturation as well as vacuum 
properties and correlate it with the star properties. 



o 
o 



> 

(N 

\o 

(N 
O 

(N 
O 

a^ 
o 



X 



I. INTRODUCTION 

Aspects of Nuclear Matter:- Nuclear matter is a 
hypothetical, infinitely large system of nucleons {N = Z) 
with coulomb interaction turned off. The matter inside 
the compact stars, the early universe, a nuclcon gas, or 
quark-gluon plasma (QGP) are some of the types of nu- 
clear matter. The only physical parameters that remains 
to characterize such a system is the bulk binding energy 
and the saturation density, letting the total number of 
nucleons go to infinity. 
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FIG. 1: Theoretical Phase diagram for nuclear matter. 



Fig. 1 shows the conventional picture of different as- 
pects and phase of nuclear matter with the variation of 
baryon density (ps) and temperature (T). At low tem- 
perature and high density, roughly of the order of ten 
times normal nuclear density, we approach the matter 
configuration which may be similar to that present in the 
core of neutron stars. The discovery of neutron stars in 
the form of pulsars has been a major stimulus to dense 
matter studies, which serves as an ideal laboratory to 
study matter at extreme conditions. They are the most 
compact objects observed till date, confined to a radius of 
nearly 10 km with the density ranging from (3 — 10) pa 
{po is the normal nuclear matter density) in the core of 
these stars. Bound by gravity, the evolution and consti- 
tution of these stars represents a beautiful amalgamation 
of all the known forces in nature. 
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Theoretically, the framework of Quantum Hadrody- 
namics P, Q present an elegant and consistent treatment 
of finite nuclei as well as infinite nuclear matter, which 
seem to provide solution to the so called "the Coester 
band^ problem 0] . A realistic nuclear equation of state 
must satisfy certain minimum criteria quantified as the 
"nuclear saturation properties", which are the physical 
constants of nature. In this context, the two most im- 
portant quantities which play vital role and are known to 
have substantial impact on the EOS are the nucleon effec- 
tive mass and the nuclear incompressibility Q , which are 
not very well determined and posses large uncertainty. 
Nucleon effective mass or the medium mass modification 
of nucleon in nuclear medium is a consequence of the 
Dirac field and forms an essential element for the suc- 
cess of relativistic phenomenology. On the other hand, 
the nuclear incompressibility derived from nuclear mea- 
surements and astrophysical observations exhibit a broad 
range of values K = (180 - 800) MeV Q. Further 
the non-relativistic and the relativistic models fails to 
agree to a common consensus. The non-relativistic cal- 
culations predict the compression moduli in the range 
K = (210 - 240)MeI/ [l], whereas, relativistic calcula- 
tions predicts it in the range (200 - 300) MeV 0, i]. 
Apart from that we are inevitably marred by the uncer- 
tainty in the determination of mass of the scalar meson 
(cr- meson) . The attractive force resulting from the scalar 
sector is responsible for the intermediate range attrac- 
tion which, along with the repulsive vector forces pro- 
vides the saturation mechanism for nuclear matter 0]. 
The estimate from the Particle Data Group quotes the 
mass of this scalar meson '/o(600)' or cr— meson in the 
range (400 - 1200) MeV [1]. A recent estimate however, 
for sigma meson mass is found to be 513 ± 32 MeV [l0| . 
Therefore in order to address these issues, one needs a 
model that has the desired attributes of the relativistic 
framework and which can be successfully applied to var- 
ious nuclear force problem both in the vicinity of po as 
well as at higher densities. 

Neutron Star Structure:- Matter at extremes opens 
up numerous possibilities, which may be present in neu- 
tron star core as shown in Fig. 2. Possible scenarios such 
as the transition from nuclear and hadronic matter, to 
exotic states involving pionic and Kaonic Bose-Einstein 
condensation, to bulk quark matter and quark matter 
in droplets, including super-conducting states, as well as 



2 




10" g/cm^ 



Strange 
star 



Nucleon 
star 



FIG. 2; Possible phases structure in neutron star core. 



strange quark matter, have been proposed 0, [HI, • Hy- 
perons can form in the neutron star cores, when the nu- 
cleon chemical potential are large enough to compensate 
the mass difference between the nucleon and the hyperon. 
Also, in dense matter, the energy of the pion (tt meson) 
and the kaon (K) is modified by inter-particle interac- 
tions and if it becomes sufficiently low, they will form 
the condensate, i.e., being bosons, they will accumulate 
the same ground state [ij] . The presence of the conden- 
sates is also known to enhance the neutrino cooling p^ . 
in addition to the nucleon direct URCA (Unrccordablc 
Cooling Agents) process [l^. It is also widely believed 
that the attractive interactions between the quarks will 
lead to pairing and color superconductivity. For three 
massless flavors, the condensation pattern that minimizes 
the free energy is known as the color-flavor-locked (CFL) 
scheme [l^ . A detailed description of these novel phases 
can be found in Ref . [13, [III . 

Observations from Neutron star:- Observations 
from neutron stars can lead to an understanding of the 
state of their interiors and the key unknowns such as the 
stars maximum mass and radius. Apart from that, the 
other key obscrvablcs are the gravitational redshift, the 
central density of the star, the moments of inertia and 
the pulsar timings. These informations can be inferred 
from the photons, ranging from radio waves to X-rays, 
and also those involving neutrinos and gravity waves. 
The largest precisely known neutron star mass is only 
1.44 Mq, but in case of binaries with white dwarf com- 
panions, the mass might be somewhere in the vicinity of 
2 Mq or even larger. The radius of the star is found to 
be sensitive to the properties of the nuclear matter near 
the saturation density (po), such as the density depen- 
dence of the nuclear symmetry energy (l9| . However, the 
radius measurements are not that precisely known unlike 
the mass observations, but the upper limits of the radius 
values are inferred from the thermal emission of cooling 
neutron stars, the gravitational redshift measurements, 
or from the crustal properties such as the pulsar glitches, 
star quakes and cooling timescales. The vital informa- 
tion from these observations can be used to constrain 
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FIG. 3: Observed mass of neutron stars in binaries. The figure 
is taken from Ref [T^j. {PSR + NS) is the binary of a pulsar 
and a neutron star and (PSR + WD) is the binary of a pulsar 
with a white dwarf companion. HMXBs' are the high-mass 
X-ray binaries and 'LMXBs' are the low-mass X-ray binaries. 
The two objects plotted as crosses, the HMXB 4U 1700 - 37 
and the LMXB 2S 0921 - 630 are blackhole candidates, but 
may be neutron stars. Error bars are la errors. 



the nuclear equation of state at high densities. 

The measured mass of 40 neutron stars in binary sys- 
tems are summarized in Fig. 3. In case of J0737 — 
3039', for the A and B system, both neutron stars are 
detected as radio pulsars, whereas in the three systems 
with the companion marked by '?' may be a white dwarf. 
In a high-mass X-ray binary (HMXB) , the companion is 
a massive star M > IOMq, whereas in a low-mass X-ray 
binary (LMXB), the companion's mass is below IMq. 
The most accurately measured masses are from timing 
observations of the radio binary pulsars such as in case 
of binary pulsar PSR 1913+16, in which the masses are 
M = 1.3867 ± O.OOO2M0 and M = 1.4414 ± O.OOO2M0, 
respectively. It is found that the mass determinations in 
binaries with white dwarf companions shows a broader 
range of neutron star masses than binary neutron star 
pulsars. The neutron star, PSR J0751-I-1807 with mass 
2.1 Mq, is about 4(t from the canonical value of 1.4 Mq. 
The system of VelaX — 1 however, has lower mass limit 
(1.6 to 1.7Af0). 

Observed radii: From the study of the measured flux 
of the nearest neutron star i?X J1856.5 — 3754, with the 
Hubble space telescope, Kaplan et.al., estimated the stars 
radius [2(|. The analysis puts the radiation radius of 
the star i?oo ~ , ^ ~ 15km corresponding to 

a / 1 — 2G AI / R 

R ~ 12 km for M ~ 1.4Af0. Recently, a few surveys 
from the non- accreting X-ray binaries, the so called the 
Quiescent stars CA:0C/132619.7- 472910.8 in ArGC5139 
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gives i?oo = 14.3 ± 2.1km (90% CL) [2i|. In case of the 
neutron star 4J71636 — 53, the compactness parameter 
'M/i?' < 0.16 implies a radius i? > 12 - 13km for the 
star with M = IAMq. 

Some other observables from compact stars arc the 
gravitational redshift measurement of the star, the mo- 
ment of inertia or the estimated central density of the 
star. Overall, there are large uncertainties in these esti- 
mates. Further stringent conditions can be imposed on 
the EOS of dense matter from the laboratory estimates of 
the heavy-ion data, such as constraints on nuclear incom- 
pressibility, symmetry energy or nuclcon effective mass. 



II. EFFECTIVE CHIRAL MODEL 

Chiral models [13, [2^ have been developed and were 
applied to nuclear matter studies. We start with an ef- 
fective chiral Lagrangian and calculate the density effects 
in nuclear matter in the mean field approach. To have 
a realistic description of the dense neutron star matter, 
we consider the effective Lagrangian of the chiral model 
generalized to include the lowest lying octet of baryons 
(ri,p, A°, S~''^) interacting through the exchange 

of the pseudo-scalar meson tt, the scalar meson a, the vec- 
tor meson ui and the iso- vector p— meson, and is given by 
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Using the meson field equations for cu and cr-meson 
along with the corresponding energy density and pressure 
expression for symmetric nuclear matter, the parameters 
of the present model arc then evaluated. Here we fix the 
saturation density (po = 0.153 fm~^) and for a desired 
value of nucleon effective mass {Y = m* /m = (0.75 — 
0.90)), we obtain the nuclear matter parameters of the 
model enlisted in Table I. For details of the procedure 
involved in obtaining the parameters, one can refer to 

We now go directly to the total energy density 'e' and 
pressure 'P' for a given baryon density in terms of the 
dimensionless variable Y = x/xq which is given as: 



e = — I k'^dkx/k'^ + m 
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The above Lagrangian represents the interaction of 
baryons b with the aforesaid mesons. We have the 
kinetic and the non-linear terms in the pseudoscalar- 
isovector pion field 'tt', the scalar field 'cr', and with 
x^ = 7?^ + cr^. Finally in the last two lines, we have the 
field strength and the mass term for the vector field 'w' 
and the iso-vector field 'p" meson. The terms in eqn. (1) 
with the subscript ' B' should be interpreted as sum over 
the states of all baryonic octets. In this paper wc shall 
be concerned only with the normal non-pion condensed 
state of matter, so wc take < tt >= 0. The interac- 
tion of the scalar and the pseudoscalar mesons with the 
vector boson generates a dynamical mass for the vector 
bosons through spontaneous breaking of the chiral sym- 
metry with scalar field getting the vacuum expectation 
value xq. Then the masses of the baryons, the scalar and 
the vector mesons, arc respectively given by 
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The terms in eqns. (3) and (4) with the subscript 'B' 
should be interpreted as sum over all the states of the 
baryonic octets. The meson field equations for the cr, uj 
and (0— mesons are then solved self-consistently at a fixed 
baryon density to obtain the respective field strengths. 
The EOS for the /3— equilibrated for the hyperon rich 
matter is obtained with the requirements of conservation 
of total baryon number and charge neutrality condition 
[2^ . Using the computed EOS for the neutron star se- 
quences, we calculate the structural properties of neutron 
stars with and without hyperon core. 



III. RESULTS 

Nuclear matter saturation is a consequence of the in- 
terplay between the attractive (scalar) and the repulsive 
(vector) forces and hence the variation in the coupling 



strength effects other related properties as well. Fig. 
mjA) reflects the same, where we have plotted the nuclear 
incompressibility for the evaluated parameter sets of the 
present model as a function of the nucleon effective mass. 
For better correlation between them, the corresponding 
ratio of the scalar and vector coupling is also indicated. 
On comparison with the incompressibility bound inferred 
from heavy ion collision experiment (HIC)[2^, we find 
that the EOS with lower nucleon effective mass is ruled 
out. The model favors EOS for which the nucleon effec- 
tive mass m* /m > 0.82. It can also be seen that the EOS 
becomes much softer with increasing ratio of C'a/C^. 
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FIG. 4: (A)- Nuclear matter incompressibility as a function 
of the nucleon effective mass for the parameters of the present 
model at po = 0.153/m~'^. Also plotted is the corresponding 
ratio of scalar to vector coupling on the opposite x-axis. (B)- 
Incompressibility as a function of obtained sigma meson mass 
for various parameter sets. The upper and the lower limit for 
incompressibility inferred from Heavy Ion Collision data (2^ 
K — (167 — 380) MeV is shown with horizontal lines. Recent 
experimental scalar meson mass limit {nicr ~ 513 ± 32) MeV 
[lOl l is depicted with vertical lines 
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FIG. 5: (A)- Sigma meson mass as a function of the pion de- 
cay constant at a fixed saturation density {po — 0.153/m~'^). 
(B)- Nucleon effective mass as a function of the pion decay 
constant. The experimental limit of the scalar meson mass 
and the pion decay constant are indicated with the horizontal 
and the vertical lines. 
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FIG. 6; Comparison of the selected EOS with the heavy ion 
collision data at high densities [2^. (A)- Comparison with 
the Symmetric Nuclear Matter data (SNM). (B)- Comparison 
with the Pure Neutron Matter data (PNM). 



Figure [UJB) shows the variation of incompressibility 
as a function of scalar meson mass obtained for various 
parameter sets. Recent experimental estimate for scalar 
meson mass rua = 513 ± 32 MeV is compared with 
the present calculation. Here we find that the EOS with 
Y ~ (0.84 — 0.86) seems to agree with the combined 
constraint from the HIC flow data and the experimental 
meson mass range. 

Figure [HI A) shows the obtained sigma mass as a func- 
tion of the vacuum value of the pion decay constant. The 
experimental bound of the pion decay constant seems to 
agree with slightly higher value of m^-, which agree with 
the upper bound of the experimental bound on nia fiot . 
FigureEI^B) shows the nucleon effective mass Y = m* /m 
as a function of the pion decay constant. The constraint 
of agree with EOS with Y = (0.82-0.84), however the 
corresponding incompressibility lies on the higher side of 
presently acceptable bounds (Q - [1]). 

Fig. OA) display the pressure as a function of baryon 
density up to nearly 6po for the selected parameters of 
the model for symmetric nuclear matter. The shaded re- 



gion corresponds to the experimental HIC data [26[ for 
symmetric nuclear matter (SNM). Among the three the- 
oretical calculations shown, the EOS with Y = 0.85 & 
0.87 agree very well with the collision data. Precisely, the 
third set (K = 231 MeV) completely agree with the flow 
data in the entire density span of 2 < pb/po < 4.6. In 
Fig. [ni^B), the case of pure neutron matter (PNM) is com- 
pared with the experimental flow data. The experimental 
flow data is categorized in terms of stiff or soft based on 
whether the density dependence of the symmetry energy 
term is strong or week [l^I- The EOS predicted by the 
present model seems to rather lie on the softer regime. 
However, the EOS with Y = 0.87, K = 231 MeV though 
satisfy the combined constraint rather well, is not consis- 
tent with the vacuum value of the pion decay constant. 

In Fig [3 we plot the results for the global proper- 
ties of the static star, obtained from the tabulated pa- 
rameters of the model for two cases. In one, the star is 
composed purely of /3— equilibrated nucleon only matter 
(n,p, e,/i) and in the other case, the star is composed 
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TABLE I: Parameter sets of the effective chiral model that satisfies the nuclear matter saturation properties such as binding 
energy per nucleon B /A — m — —16.3 AleV, nucleon effective mass Y = m* /m = (0.75 — 0.90) and the asymmetry energy 
coefficient is J « 32 MeV at saturation density po = 0.153/m~''. The nucleon, the vector mesons {u> & po) masses are taken to 
be 939 MeV, 783 MeV and 770 MeV respectively and — {g„/m^Y , = {g^/m^Y s^nd Cp = {gp/'nipy' are the corresponding 
coupling constants. B = b/m^ and C = c/rn^ are the higher order constants in the scalar field. Also given is the scalar meson 
mass 'mo^', the pion decay constant '/^' and the nuclear matter incompressibility [K) at po- The maximum mass and radius 
of neutron star composed of nucleon only matter {Mm, Rn) and hyperon rich matter {Mr, Rh) for static case are also 
tabulated. The last column shows the central density of the star composed of nucleon [en) and hyperon (en)- 



set c„ Cuj Cp B C rria Y f^r K Mn{Mh) Rn (Rh) £n (en) 

ifm^) ifm^) (/m^) (/m^) (/m^) (MeV) (MeV) (MeV) (M/Mq) (km) (lO^^gem^) 

1 5.916 3.207 5.060 1.411 1.328 691.379 0.75 110.185 1098 2.49 (2.04) 14.2 (20.8) 1.26 (0.59) 

2 6.047 3.126 5.087 0.822 0.022 675.166 0.76 111.601 916 2.47 (2.02) 14.2 (20.7) 1.26 (0.62) 

3 6.086 3.031 5.107 0.485 0.174 662.642 0.77 113.346 809 2.44 (1.98) 14.0 (20.5) 1.30 (0.62) 

4 6.005 2.933 5.131 0.582 2.650 656.183 0.78 115.238 737 2.41 (1.95) 13.9 (20.4) 1.35 (0.62) 

5 6.172 2.825 5.155 -0.261 0.606 635.287 0.79 117.403 638 2.38 (1.92) 13.7 (20.2) 1.40 (0.65) 

6 6.223 2.709 5.178 -0.711 0.748 619.585 0.80 119.890 560 2.34 (1.89) 13.6 (20.0) 1.44 (0.65) 

7 6.325 2.585 5.200 -1.381 0.089 600.270 0.81 122.740 491 2.30 (1.85) 13.4 (19.8) 1.49 (0.68) 

8 6.405 2.451 5.222 -1.990 0.030 580.876 0.82 126.039 440 2.26 (1.81) 13.2 (19.5) 1.54 (0.68) 

9 6.474 2.323 5.242 -2.533 0.300 562.500 0.83 129.465 391 2.22 (1.77) 13.0 (19.3) 1.60 (0.71) 

10 6.598 2.159 5.265 -3.340 0.445 536.838 0.84 134.378 344 2.16 (1.72) 12.8 (19.0) 1.65 (0.71) 

11 6.772 1.995 5.285 -4.274 0.292 509.644 0.85 139.710 303 2.10 (1.66) 12.6 (18.7) 1.76 (0.75) 

12 7.022 1.823 5.305 -5.414 0.039 478.498 0.86 146.131 265 2.04 (1.62) 12.3 (18.3) 1.82 (0.78) 

13 7.325 1.642 5.324 -6.586 0.571 444.614 0.87 153.984 231 1.97 (1.56) 12.0 (17.9) 1.94 (0.82) 

14 7.865 1.451 5.343 -8.315 0.502 403.303 0.88 163.824 199 1.89 (1.51) 11.8 (17.5) 2.08 (0.86) 

15 8.792 1.249 5.362 -10.766 0.354 353.960 0.89 176.552 168 1.80 (1.44) 11.4 (16.9) 2.30 (0.94) 

16 7.942 1.041 5.388 -6.908 15.197 339.910 0.90 193.437 163 1.71 (1.37) 11.0 (16.3) 2.54 (1.03) 



e (g cm ^ ) 

S.Oxlo"* l.Oxlo'' 1.5xl0'^ 2.0xlo'' 2.5xlo'" 
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FIG. 7: ( Upper panel:-) The mass of the star obtained as a 
function of the central density of the star. {Lower panel:- 
) The mass of the star as a function of the radius of the 
star. The two cases considered here are the star composed of 
nucleons only (n, p, e, /x) shown with the filled circles and the 
one composed with hyperons, shown with the triangles. 



the vector counterpart to yield the binding of A as sug- 
gested by hypernuclei experiments. In the upper panel 
of the figure, the maximum mass of the star is plotted 
as a function of the central density of the star. Here it 
can be seen that the star with nucleon matter is much 
massive than those composed of hyperons. However in 
case of star composed of nucleons, the central density 
ranges from (4.5 — 9.5) x po, but with hyperons, the cen- 
tral density of the star falls in the range « (2 — 4)x po 
and the resulting radius increases for the later case. For 
all the cases given in Table I, we find nearly (46 — 48) 
% increase in the star radius, when we move from nucle- 
onic star to hyperon rich star matter. The corresponding 
decrease in maximum mass is « (18 — 20)%. This fea- 
ture is very much evident from the lower panel of Fig. [3 
where we have plotted the maximum mass of the star ob- 
tained as a function of the star radius. In the transition 
of nucleon only matter to hyperon rich neutron star mat- 
ter, the compactness ratio (M/R) of the star falls from 
(0.15-0.17) to (0.09-0.10). 



IV. SUMMARY 



of matter which is Hyperon rich. To study hyperon rich 
matter, we fixed the hyperon couplings for all the pa- 
rameter sets, so as to yield the binding energy of A° at 
Po = —30MeV. This is obtained by fixing the scalar cou- 
pling Xa- = XaN /xaH = 0.70 and varying the strength of 



The effective chiral model provides a natural frame- 
work to interlink the standard state properties of nu- 
clear matter with the vacuum correlations, such as the 
pion decay constant. We find that the pion decay con- 
stant is experimentally well known quantity in compar- 
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ison to other derived quantities such as the nuclear in- 
compressibiUty and cr— meson mass, that can put strin- 
gent constraint on the model parameters. Experimen- 
tally determined effective mass from scattering of neu- 
tron over Pb nuclei [sH seems to go well with the present 
model. Both of them favor higher value for nucleon ef- 
fective mass. On a comparative analysis of the resulting 
EOS with that of the HIC data for symmetric nuclear 
matter as well as pure neutron matter, parameter set 
with Y = 0.85; K « 300 MeV seems to be the ideal pa- 
rameterization of the present model. The resulting scalar 
meson mass mg- « 510MeV, is also consistent with the 
experimentally observed masses [l^, [HI . We then applied 
the model to study the global properties of neutron stars 



with and without hyperons. We find that the star com- 
posed of hyperons has a lower value of central density and 
results in larger radius than the star without hyperons. 
This in turn results in the decrease in the compactness 
parameter of the star. Although the resulting higher cen- 
tral density of the nucleonic star seems to be unphysical 
on account that at those higher densities, it is very un- 
likely that nuclcons would retain its identity. Further, 
from the astrophysical point of view (observed neutron 
star masses), none of the parameters can be ruled out, 
however it shall be interesting to study the rotational at- 
tributes of the star with hyperons/ mixed/ quark phase 
in order to constrain the parameters of the model. Work 
is in progress in this direction (s^l • 
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